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In this paper, spectral analysis of fractional Sturm Liouville problem defined on (0, 1], having the singularity of 
type — - + at zero and research the fundamental properties of the eigenfunctions and eigenvalues for the 

operator. We show that the eigenvalues and eigenfunctions of the problem are real and orthogonal, respectively. 
Furthermore, we give some important theorems and lemmas for fractional hydrogen atom equation. 
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(— | • 1. Introduction 

Sturm-Liouville problem was first developed in a num- 
^nber of papers that were published by these authors in 1836 
+^and 1837. (Charles-Francois Sturm (1803-1855), Professor 

S of Mechanics at the Sorbonne, had been interested, since 
about 1833, in the problem of heat flow in bars, so was well 

aware of eigenvalue- type problems. He worked closely with 
^Hhis friend Joseph Liouville (1809-1882), Professor of Mathe- 
matics at the College de France, on the general properties of 
^~^s,econd order differential equations. Liouville also made many 
qq contributions to the general field of analysis.) [1]. Mathe- 
^Njmatical techniques could be developed into a more suitable 
and significant course by presenting them within the more 
general Sturm-Liouville theory in L2. The Sturm-Liouville 
CO Problems are important in many areas of science, engineer- 
^Hing and mathematics. It is known that the spectral char- 
acteristics are spectra, spectral functions, scattering data, 
norming constants, etc. According to the theory a linear 
second-order differential operator which is self-adjoint has an 
^orthogonal sequence of eigenfunctions in L^. Spectral prop- 
erties of Sturm-Liouville operators are often derived, directly 
or indirectly, as a consequence of an established link between 
large distance asymptotic behavior of solutions of the associ- 
ated differential equation and spectral properties of the cor- 
responding differential operator. Sturm-Liouville problems 
are divided into regular and singular type. Differential equa- 
tions such as Bessel, hydrogen atom, Hermitte, Jakobi, and 
Legendre equations can be transformed into Sturm Liouville 
equations. There are many studies on these issues, [1-5] . We 
also discuss the radial part of Schrodinger's equation for the 
hydrogen atom. 

Fractional calculus is 'the theory of derivatives and in- 
tegrals of any arbitrary real or complex order, which unify 



X 



and generalize the notions of integer-order differentiation and 
n— fold integration' [6-13]. In recent years, the concept of 
fractional calculus, originated from Leibniz, has achieved in- 
creasing interest during last two decades. In particular, on 
the last decade have scientific papers concerning fractional 
quantum mechanics. It has been proved that many sys- 
tems in different fields of science and engineering can be 
modeled more accurately using fractional derivatives [14,15] . 
We note that ordinary derivatives in a traditional Sturm- 
Liouville problem are replaced with fractional derivatives and 
the resulting problems are solved using some numerical meth- 
ods [16-20] . Furthermore, Klimek and Agrawal [21] define a 
fractional Sturm-Liouville operator, introduce a regular frac- 
tional Sturm-Liouville probelem and investigate the proper- 
ties of the eigenfunctions and the eigenvalues of the operator. 
In this paper, our purpose is to introduce regular fractional 
singular Sturm-Liouville problem and prove spectral proper- 
ties of spectral data for the operator. 

Let's give the main problem and necessary datas, as fol- 
lows: 

2. Preliminaries 

Firstly, we consider singular Sturm-Liouville equation 



d 2 R 
dr 2 



a dR 
r dr 



£(£ + !) 



R+[E 



R = 



(0 < r < 00) . 



In quantum mechanics the study of the energy levels of the 
hydrogen atom leads to this equation [5] . Where R is the dis- 
tance of the mass center to the origin, £ is positive integer, 
a is real number E is energy constant and r is the distance 
between the nucleus with the electron. 

Substitution R — y/r reduces above equation to the form 
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where z,ai,bj G C,cti,0j G R(i = 1, ...,p;j = 1, ...q) . This 
^ ^ + { E + a ^ ^ 1 — general Wright function was investigated by Fox who pre- 

dr 2 \ r r 2 J ^ sented its asymptotic expansion for large values of the argu- 

ment z under the condition 



In most cases we took the reference potential to be 

a £(£+1) 



q p 



- J2 ai > L 



If these conditions are satisfied, the series in (7) is convergent 
when V (r) is defined by last differential equation, it contains for any z G C. 

a centripetal and Colomb part, the usual singularities of the Theorem 4. [14] Let ai,bj G C, ai,/3j G R (i = 1, ...,p) , 
nuclear problem. (.7 = 1, ■••«) and let 

q p 

Definition 1. Let Re (a) > 0. The left-sided and re- A = — y^ctj, 

spectively right-sided Ricmann-Liouville integrals of order a j=i i=i 

are given by the formulas p „ 9 fl 

(/« + /) (r) = j^j (r - 1 / W r > a, (1) ^ = £ fc . _ £ ^ + p-. 



I) If A > —1, then the series in (7) is absolutely conver- 
gent for all z G C. 

II) If A = —1, then the series in (7) is absolutely conver- 
gent for \z\ < S and for \z\ = S and 5ft (p) > 1/2. 

where L denotes the gamma function. Definition 5. [22] Let (X, d) be a non-empty complete 

Definition 2. Let Re (a) > 0. The left-sided and re- metric space. Let T : X -> X be a contraction mapping on 
spectively right-sided Riemann-Liouville derivatives of order r, i.e.. there is a nonnegative real number g < 1 such that 
a are defined as d (T (r) ^ T < ^ (/ , y) 

(D° + /) (r) = D +" f) (r) r > a, (3) f° r au r , 2/ in A. Then the map T admits one and only one 

fixed-point r* in X. (this means T(r*) = r*). Furthermore, 

(D a f)() D ( I 1 ^" f \ ( b (4) ^ S ^ xc< ^ P om * can ^ e f° un< i as follows start with an arbi- 

( b-J ) v) ~ \ b- J J \ r ) r < > w trary element r in r and define an iterative sequence by r n 

, „ , . i, ,i in i • , , • i , r~i = T(r r .--[ ) for n = 1, 2, 3, ... This sequence converges, and its 

analogous formulas yield the left- and right-sided Caputo ... . ' _ n . „ . . ,. , ., , , r 

j ... r j ' limit is r . I he following inequality describes the speed of 

derivatives of order a: ° ^ J r 

convergence: 

(r) = (i^Df) (r) r > a, (5) d(rVn) < ^(n.ro) . 

1 — q 

Equivalently, 

(^V/)W = (/;-»(- C )/)W r< b . (6) ^rjs^rj 

1 — <J 

Definition 3. [14] The general function (z) is de- anc j 
ed for z G C, complex ai,bj G 
. (i = 1, ...,p;j = 1, ...<?) by the series 



fined for z G C, complex aj.bj G C, and real q,,^ t d (r*, r„ + i) < erd (r*, r„) . 



j>*9 (z) =„ * g 



(ai 


«i)i. P 




. ( & i 




z 



Any such value of a is called a Lipschitz constant for T, 
and the smallest one is sometimes called "the best Lipschitz 
constant" of T. 

Property 6. The fractional differential operators de- 
Y[ r (cn + aik) fined in (3 — 6) satisfy the following identities: 

ysi — — (7) b b 

k=o n r (bj + frk) jf (r) Dl_g (r) dr = J g (rf D° + / (r) dr-f (r) /^.g (r) |* 



P 

OO 



2 



or 





j f{r)Dl_g{rf D a a+ k{r)dr 

a 

b 

= jg{r) C Dl + f{r) C Dl + k{r)dr (8) 

a 

- !{r)ll_ a g{rfDl + k{r) \ \, 

b 

J f{r)Dl + g{r)dr 

a 

b 

= J g (rf Dl_f (r) dr + f (r) I^g (r) \ b a . (9) 



Property 7. Assume a £ (0, 1), /3 > a and / £ C [a, b] . 
Then the following relations 



Da, + Ia,+f(r) = f(r) 
Dl + li + f{r) = lt + a f{r) 



lL a f{r) 



^,+^,+/W=/(r) 
c D? i _I£_f(r)=f(r) 



(10) 



hold for any r e [a, 6]. Furthermore, the integral operators 
defined in (1,2) satisfy the following semi-group properties, 



J. 



a+/3 
o,+ ' 



1 b,- 1 b,- 



r«+/3 



Now, let's take up a fractional Sturm-Liouville problem 
for hydrogen atom equation 

3. A Singular Fractional Sturm-Liouville Problem 
for Hydrogen Atom Equation 

Let's denote a singular fractional Sturm-Liouville prob- 
lem for hydrogen atom equation with the differential part 
containing the left-and right-sided derivatives. Let's use the 
form of the integration by parts formulas (8,9) for this new 
approximation. Main properties of eigenfunctions and eigen- 
values in the theory of classical Sturm-Liouville problems are 
related to the integration by parts formula for the first order 
derivatives. In the corresponding fractional version we note 
that both left and right derivatives appear and the essential 
pairs are the left Riemann-Liouville derivative with the right 
Caputo derivative and the right Riemann-Liouville derivative 
with the left Caputo one. 



Definition 8. For a £ (0, 1) , Fractional hydrogen atom 
operator is written as 



- — + qo(r)) (11) 



consider the fractional hydrogen atom equation 



£ a y\ (r) + Xw a (r) y x (r) = 



(12) 



where p (r) ^ 0, w a (r) > Vr £ (0, 1] and p, q are real val- 
ued continuous functions in interval (0,1]. The boundary 
conditions for the operator C a are the following: 



yx (o) = o, 



(13) 



diy (1) + d 2 l\;_ a p (if Dl +y (1) = 0, (14) 



where d\ + d\ ^ 0. The fractional boundary-value problem 
(12)- (14) is fractional Sturm-Liouville problem for hydrogen 
atom operator. 

Theorem 9. Fractional hydrogen operator is self-adjoint 
on (0, 1] . 

Proof. Let us consider the following equation 



(C a ip,(j>) = J C a ip(f>dr 

o o 
+ , 2 t 1(1 + 1) 



Dl_p(r) c Dl +V (r) 



+ 9o (r) V (r) 



dr 



J<t>(r)Dl_p(rf D% + cp(r)dr 



2 1(1 + 1) . . . 
- + 2 + qo (r) )f(r)(f> (r) dr. 



Considering property (8) to the first integral in the last 
equation, we obtain the identity 
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1 

(£ a <p, ^)=Jp (rf D^ + cb (rf Dg >+ <p (r) dr 



-<P(r)ll_y(rf Dl + <p{r) 



2 1(1 + 1) 

— + 2 + 9o (r) ) y (r) (r) dr 



1 

= Jp(rf D% !+ <j>(rf D% !+i p(r)dr 



+ % (1)0(1) + % (0)0(0) 

«2 C 2 



2 + 



9o (»") ) y (r) (r) dr. (15) 



Similarly, we obtain 



and its complex conjugate y 

C a y(r) + Xw a (r)y(r) =0 (18) 

JA (0) = (19) 

d W (1) + da/^p (1) C Dl + y (1) = (20) 

C a y(r)+Xw a (r)y(r) = (21) 

W(0) = (22) 

diy (1) + d^ 1 ""? (1) C Dl + y (1) = (23) 

where d\ + d\ ^ 0. We multiply equation (18) by function y 
and (21) by function y respectively and subtract: 

(A - A) w a (r) y (r) y (r) = y (r) £ Q y (r) - y (r) £ Q y (r) . 

Now, we integrate over interval (0, 1] and applying relation 
(17) we note that the right-hand side of the integrated equal- 
ity contains only boundary terms: 



(<P, ^4>) = jp (rf (rf D£ + (r) dr 

o 

+ ^(1)0(1) + % (0)0(0) 



1 

(A - A) J w a (r) y (r) y (r) dr 



.2 + /([+l) +9o(r))()(/ . ) ^ (/ . )r/r (;|()) 



The right hand sides of the equations (15) and (16) are equal 
hence we may see that the left sides are equal that is 

(C a Lp,(f>) = (ip,C a (/)) . 

Therefore, C a = C* a . The proof is completed. 

Theorem 10. The eigenvalues of fractional hydrogen 
atom operator (12 — 14) are real. 

Proof: Let us observe that following relation results from 
Property (8) 

l l 

jf (r) C a g (r) dr = jp (rf D« + f (rf D« + g (r) dr 



y (r) C a y (r) dr - y (r) C a y (r) dr 



1 

Jy (r) 



Dl_p(rf Dl + y(r) 



+ — + 



2 1(1 + 1) 



+ qa (r) y (r) 



dr 



Dl_p(rfDl + (r)y(r) 
2 1(1 + 1) \ 

- + r 2 + 9" ( r ) ) v ( r > 

± 

Jp(rf Dl + y(rf D« + y(r)dr 



-f(r)l\-«p(rf Dl + g(r) 



-y(r)ll-_ a p(rf ' (r) 



+ 



l + l l+ll + qo{r) ) fj{r] f [r](ll , , 17) 



+ 



/(- 



2 1(1 + 1) 



+ 



+ Qo (r) )y(r)y (r) dr 



Suppose that A is the eigenvalue for (12 — 14) correspond- 
ing to eigenfunction y the following equalities are satisfy y 



i 

Jp(rf Dl + y(rf Dl + y(r)dr 



4 



or 



+y(r)Il-_ a p(r) C D% !+ y(r) 



o 



Integrating over interval (0, 1] and applying relation (17) we 
note that the right-hand side of the integrated equality con- 
tains only boundary terms: 



2 1(1 + 1) , , , 
- H 2 \-qo(r) ) y(r)y(r)dr 



(A -A) jw a (r)\y(r)\ 2 dr 
o 

= -y{r)l{;_ a p{rf D° + y(r) 
+ y(r)ll-_ a p(rf D« + y(r) 
+y{r)l\;_ a P {rf Dl + y{r) 
-y(r)ll-_ a p(rf D« + y(r) 



by virtue of the boundary conditions (19) , (20) , (22) , (23) , 
we find 

l 



(Ai - A 2 ) J w a (r) y Xl (r) y\ 2 (r) dr 
o 

l l 

= J VXi ( r ) £ a y\ 2 (r) dr - j yx 2 (r) C a y Xl (r) dr 



1 



Dl_p{rf Dl + y X2 (r) 



+ 



2 1(1 + 1) 



+ 



+ qo (r) y\ 2 (r) 



dr 



- J VX2 ( r 



r) Dl_p{rf D« + (r)y Xl (r) 



( a-a)/ 



w a(r)\y(r)\ dr = 



+ 



2 1(1 + 1) 



+ 



+ Qo (r) y Xl (r) 



dr 



and because y is a non-trivial solution and w a (r) > 0, it 
easily seen that A = A. 

Theorem 11. The eigcnfunctions corresponding to 
distinct eigenvalues of fractional hydrogen atom operator 
(12 — 14) are orthogonal weight function w a on (0,1] that 
is 



/ 



w a (r) y Xl (r) y\ 2 (r) dr = Ai ^ A 2 . 



Proof. We have by assumptions fractional Sturm- 
Liouville for hydrogen atom operator fulfilled by two dif- 
ferent eigenvalues (Ai,A 2 ) and the respective eigenfunctions 
(yx!,yx 2 ) , m tnat case 

C a y Xl (r) + X lWa (r) y Xl (r) = (24) 

Vx 1 (r) = (25) 

diy Xl (1) + d 2 l\:_ a P (if Dl +VXl (1) = (26) 

C a y X2 (r) + X 2 w a (r) y X2 (r) = (27) 

yx 2 (r) = (28) 

diy X2 (1) + d 2 l\:_ a v (if Dl +VX2 (1) = (29) 

we multiply equation (24) by function y X2 and (27) by func- 
tion y Xl respectively and subtract: 

(Ai - A 2 ) w a (r)y Xl y X2 = y Xl C a y X2 - y x . 2 C a y Xl . 



Using the boundary conditions (25) , (26) , (28) , (29) , we find 
that 

l 

(Ai - A 2 ) jw a (r) y Xl (r) y X2 (r) dr = 0, 
o 

where Ai ^ A 2 . 

Remark. Let us now give certain auxiliary functions. 
Because we use the functions, the first of them is as follows 

xct-l 



lo j 1 ^!, =(i-or i (r-or 



r» 



•1 *2 



(1.1) 

(a,-l) (a + 1,1) 



r - 
1-0 



where i^ , 2 is the Fox- Wright function [14]. 

1*2 

r (ai + ctik) 



(ai,ai) 

(Mi) (b 2 ,/3 2 ) 



E 

fc=0 



r(6i +/3 1 k)T(b 2 +p 2 k) k\ 



the properties of the function are determined by the param- 
eters 

A = /3i +/? 2 -ai = -1 



1-2 „ 1 

H = bi + b 2 - ai H — = 2a - -, 



(30) 



considering Theorem 4, we note that this function is con- 
tinuous in (0,1] when order a > 1/2, i.e. \i > 1/2. It is 
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discontinuous at end r = 1 for < a < 1/2. The explicitly 
calculated function allows to estimate the second component 
of stationary function </> of the differential part of diffusion 
operator 

Dl_p{rf Dl + <p (r)=0 
which looks as follows 

00 (r) = a + 6J °+ ( r Vw ^ = & + ^ °> r ) • ( 31 ) 

' 1 (a) p (r) 
The next function is the following integral 

^ (r) = /oV"- 1 ^ 7 o >+ ^^ = (! " °)" ( r - °)" 



Xi ^ 



(1.1) 

(a + 1,-1) (a + 1,1) 



r -0 



1-0 



(32) 



Again, using the Theorem 4 and calculating parameters ac- 
cording to (30) 



A = -l, 5 = 1, ^ = 2a+i. 



we conclude 



a > 



and the obtained Fox- Wright function (32) is continuous 
interval (0, 1] for any positive order a. 
Lemma 12. Let a > 1/2 and denote 



m 



Yx(y) 



-- + 2 + 9o (r) 2/a (r) + A™ Q y A (r) 



A = rf 2 + di^(a,0,l). 



Assume A + 1 0. Then, on the C [a, 6] -space, the regular 
(12 — 14) is equivalent to the integral equation 



Vx {r) = -I% +J faI? t _Y x (y) 
+Mr) (lS,+ ttI?,-Y\(vj) lr=i 



(33) 



where the coefficient A (r) is 

A(r) = |-V(a,0,r) 

and functions V> is defined in (31) . 

Proof. By aid of composition rules, equation (12) can 
be rewriten as follows: 



Dl_p(r) c Dg >+ 



1 



= 0. 



The last equality refers that on the C [a, b] -space the function 
in brackets is a stationary function of fractional Sturm Liou- 
ville problem for hydrogen atom operator. Df \_p(r) C Dq + 
which according to (31) can be found as 



Si + 42J ,+ - 



Ci+6^(a,0,r) 



T(a)p(r) 
equation (12) in the form of 

yx (r) + IZ,+^I?,-Yx (y) = 6 + & V («, 0, r) , (34) 

to end the proof we should connect coefficients £j> j = 1,2 
values dj, j = 1,2 determining the boundary conditions 
(13 — 14). Let us note that the following formula results 
from composition rules (9) and equation (34) 

I\;_ a p (rf Dl +yx (r) = -ll_Y x (y) + 

for continuous function y x , we obtain the following values a 
the ends 



Ifpirr D% t+ y x (r) 



r=0 



o 



t>W flo,+» (r) 



r=l 



respectively for y x , using (34) we find 

yx (0) - 0o (0) = a 



yA (l) = ^ (l)-/« + — /?_y A ( y ) 



6 +6^ («, 0,1)-^+^/^ ( y ) 



(35) 
(36) 

(37) 
(38) 



The following set of linear equations for coefficients a results 
from (35 - 39) 

a = o 

d 1 £ 1 +b{d 2 + d 1 ip{a,0,l)) = d 1 F (39) 
where F = Ig ^If _Y x {y) . 

f\ / r—l 

Since A ^ 0, the solution for coefficients a = 1, 2) is 
unique: 

a = o 

a = diF/L 

substituting the above solution into (34) we recover the 
equivalent integral equation (33) . 

Furthermore we give notation such as 



m, D = mm 



\P Ml . 



r€(0,l] 

A=\\A(r)\\, M v = y(r)\\, 
where ||.|| denotes the supremum norm on space C [0, 1] .( 



G 



Theorem 13. Let a > 1/2. Suppose that A ^ 0. 
Then unique continuous eigenfunction y x of singular frac- 
tional Sturm-Liouville problem for hydrogen atom operator 
(12 — 14) corresponding to each eigenvalue obeying 



2 1(1 + 1) 



+ q (r) + Xw a 



< 



M v + Bip (1) 



(40) 



on space (C [0, 1] , ||.||) 

\\Tg-Tp\\<\\g-p\\L, L G (0, 1] . 

Hence, a unique fixed point enounced as y x G C (0, 1] exists 
that solves equation (12,33) and satisfies the boundary con- 
ditions (13, 14), provided (40) is applied. In that case such 
eigenvalues are simple. The proof is completed. 



exists and such eigenvalue is simple. 

Proof. We have to say that equation (33) can be inter- 
preted as a fixed point condition on function space C [0, 1] , 

V\ = Ty x 

where the mapping on the right-hand side is for any contin- 
uous function g G C [0, 1] defined as 



Tg = -^W) 11 -^ {y)+A {r) ^ I "- + ^ It - Yx {y) 



~p(r)~ 



Conclusion 

In this article, we have extended the scope of some proper- 
ties of fractional Sturm-Liouvillc problem for hydrogen atom 
operator. We pointed that its eigenvalues related to the frac- 
tional Sturm-Liouville problem for hydrogen atom operator, 
with the certain boundary conditions are real and its eigen- 
functions corresponding to distinct eigenvalues are ortohogo- 
nal. Furthermore, we showed that the operator is self adjoint, 
l Our new results are original and important for the fractional 
Sturm-Liouville theory. 



The following inequality will be useful in further estimations: 

2 1(1 + 1) 



\\Y x (g)-Y x (p)\\ = 
+ \w a g) - 

(g-p) 



i(i 



r 
1) 



Qo (r) g 



2 

— + 
r 



1(1+1) 



+ Qo (r) I p + Aw Q P 



+ q (r) + Xw, 



< \\g-p\\ 



i(i + i) 



+ 9o (r) 



By performing necassary operations for distance between 
images Tg and Tp for a pair of arbitrary continuous functions 
g,peC[0,l]. 



\\Tg-T P \\ 
+A (r) 



1 



1 



°- + p(r) 
1 



Ii.-Yx(g) 



0,+ p{r y 

1 



II_Y X (p) +A(r)[ I« + —I«_ Yx (p) 



+A(r) 

<\\g-p 
< h-p\\L 

where L = 



AY x (g)-Y x (p)) 
Il_ (Y X (g) Y X (p)) 
9o (r) 



p(r) 

2 1(1 + 1) 



XlVr 



M, 



m„ 



m„ 



By means of (40) , we state that mapping T is a contraction 
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